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METHOD FOR MEASURING THE ABSORPTION COEFFICIENT 
AND THE REDUCED SCATTERING COEFFICIENT 
OF A MULTIPLE SCATTERING MEDIUM 
BACKGROUND OF THE INVENTION 
5 1. Field of the Invention 

The invention relates to a method for measuring the 
absorption coefficient and the reduced scattering 
coefficient of a medium, more particularly to a method 
for measuring in real time the absorption coefficient 
10 and reduced scattering coefficient of a multiple 
scattering medium. 

2. Description of the Related Art 

Human tissue is a highly scattering medium, and has 
the characteristics of having a scattering coefficient 

15 much greater than the absorption coefficient thereof. 

A conventional method of obtaining an image in a 
scattering medium involves measuring of absorption 
coefficients so as to obtain an image with contrast. 
However, when it is desired to obtain the image of an 

20 object in a multiple scattering medium, since the light 
waves are highly scattered, the resultant image is 
blurred, and the resolution of the image is also reduced 
considerably. Therefore, if the scattering effect of 
light in a relatively high scattering medium can be 

25 reduced, or if only slightly scattered snake photons 

and ballistic photons are selected, together with 
enhancement of the sensitivity to absorption 
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coefficients, the resolution of an image in the medium 
can be enhanced. However, this is not suitable for 
imaging objects in multiple scattering media. 

At present/ methods for imaging in a multiple 
5 scattering medium mainly include time-domain and 

frequency-domain imaging techniques- The concept of 
diffused photon density wave (DPDW) is proposed in the 
frequency-domain technique. DPDW satisfies the 
diffusion equation, and can be relied upon to obtain 

10 definite amounts of the absorption and scattering 
coefficients of the test object, thereby permitting 
recovery of the image of an object in the scattering 
medium. The frequency-domain technique is currently 
more suitable for imaging in a multiple scattering medium, 

15 but has the drawback that the spatial resolution is not 
high- Thus, how to enhance imaging resolution in 
applications that involve a multiple scattering medium 
is currently an important topic in the industry. 
SUMMARY OF THE INVENTION 

20 A two-frequency polarized laser , e.g., aZeemanHe-Ne 

laser, can produce a two- frequency mutually correlated 
and mutually orthogonal linear polarized photon pair 
(PPP) laser beam for propagation through a highly 
concentrated scattering medium to produce a diffused 

25 polarized photon pair density wave (DPPDW) . 

Particularly, by means of a polarizer, the mutually 
correlated and mutually orthogonal linear polarized 
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photon pair can be effectively converted into a mutually 
correlated and mutual ly paral lei polari zed photon pair . 
The present invention primarily contemplates the 
propagation of a mutually correlated and mutually 
5 parallel two- frequency polarized photon pair in a 
multiple scattering medium to form a diffused polarized 
photon pair density wave while satisfying the diffusion 
equation. The same result can be achieved if a 
two- frequency circular polarized photon pair is 
10 converted into a mutually parallel polarized photon pair 
by the use of a polarizer. By measuring directly the 
amplitude and phase of the DPPDW, the absorption 

coefficient //^ and the reduced scattering coefficient 
H\ of the multiple scattering medium can be measured 

15 in real time, in which //^ = /i^(l-g^) , where //^ is the 

scattering coefficient, and g is the scattering 
anisotropic parameter of the multiple scattering medium 

The three-dimensional (3-D) distributions of ^\ and 

//^ are obtained by finding the solution to the diffusion 

20 equation so as to obtain optical characteristics of and 
images in the multiple scattering medium. Further, as 
oxyhemoglobin (Hb02) and deoxyhemoglobin (Hb) have 
different absorption coefficients with respect to the 
same wavelength, by using a two- frequency laser source 

25 (preferably, two two- frequency lasers) having 
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different center wavelengths, the hemoglobin 
saturation (Sa02) in blood can be measured in real time, 
and images of hemoglobin saturation distribution can 
also be obtained. 
5 This invention provides a method using a 

two- frequency, two-polarization laser beam for 
propagation in a highly concentrated scattering medium 
so as to form DPPDW. Through optical heterodyne 
interference techniques and coherence of the polarized 
10 photon pair, the magnitude of the phase and amplitude 
of the DPPDW can be obtained from the heterodyne 
interference signals, and the absorption coefficient 

//^ and the reduced scattering coefficient //^ of the 

scattering medium can be calculated from the phase delay 
15 and amplitude attenuation that change with the distance 
between a laser source and a photo detector. 

The object of this invention is to provide a method 
for measuring the absorption coefficient and the reduced 
scattering coefficient of a multiple scattering medium, 
20 in which, due to coherence of a polarized photon pair, 
and processing with the use of a heterodyne polarized 
interferometer, scattering of a light beam by the medium 
can be reduced, thereby enhancing resolution. 

Another object of this invention is to provide a 
25 spatial imaging method, in which, by changing measuring 

positions and recording the measurement results at each 
of the measuring positions, a complete image of the 
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multiple scattering medium can be constructed. 

According to the present invention, the method for 
measuring the absorption coefficient and the reduced 
scattering coefficient of a multiple scattering medium 
5 includes : outputting a coherent light beam, the coherent 

beam including linear polarized P and S wave components 
having mutually orthogonal polarizations and 

frequencies co p and G) s, respectively; splitting the 

coherent light beam into a signal beam and a reference 

10 beam, the signal beam and the reference beam including 
the P wave and S wave components, respectively; 
projecting the signal beam into the multiple scattering 
medium; detecting and converting an optical 
interference signal of the reference beam and an optical 

15 interference signal of the signal beam that penetrates 

the multiple scattering medium and a polarizer, 
respectively, into heterodyne interference electrical 
signals; comparing the two heterodyne interference 
electrical signals to obtain amplitude attenuation and 

20 phase delay of the signal beam that penetrated the 
multiple scattering medium; and inferring the reduced 
scattering coefficient and the absorption coefficient 
of the multiple scattering medium at a position where 
the multiple scattering medium is penetrated with 

25 reference to the amplitude attenuation and phase delay 
thus obtained. 

BRIEF DESCRIPTION OF THE DRAWINGS 



6 



408 338427 



Other features and advantages of the present 
invention will become apparent in the following detailed 
description of the preferred embodiment with reference 
to the accompanying drawings, of which: 

Figure 1 is a block diagram showing the construction 
of an apparatus for implementing a preferred embodiment 
of a method according to the invention; 

Figure 2 is ' a schematic view illustrating the 
principle of generating a photon pair with parallel 
polarization directions in the apparatus of Figure 1; 

Figure 3 is a schematic view of an imaging device 
employed in the preferred embodiment; 

Figure 4 is a block diagram of another apparatus for 
implementing the preferred embodiment; 

Figure 5 illustrates how P wave and S wave components 

X 

are adjusted using a — wave plate in the apparatus of 

Figure 4 ; and 

Figures 6A and 6B respectively illustrate the 
variation in amplitude attenuation and phase delay of 
DPPDW obtained in an experiment conducted by the 
appl icant . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Referring to Figure 1, a two- frequency , 
two-polarization laser 100, such as a Zeeman laser, 
outputs a coherent light beam, which includes mutually 
correlated and mutually orthogonal linear polarized 
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waves that include a P wave and an S wave, the temporal 
frequencies of which being cOp and cOg, respectively. The 
coherent light beam is passed through a polarizer 200 

having an azimuth angle 0 with respect to an axis (X) . 

As shown in Figure 2, since the angle of polarization 
direction of the polarizer 200 relative to the P wave 

is 6 , the component of the P wave projected in the 
polarization direction is = cos0cxp[ia}pt) , whereas the 

projected component of the S wave is = A^sin0cxp{ia}^t) , 

thereby generating a photon pair (Lp, Lg) with mutually 
parallel polarization directions, wherein Ap and As are 
the amplitudes of the P wave and S wave components, 
respectively. Thus, a two- frequency mutually 

correlated and mutually parallel linear polarized 
photon pair beam is generated accordingly. 

Abeam splitter 210 splits the linear polar i zed phot on 
pair beam into a signal beam Ig and a reference beam Ir. 
The signal beam Ig is focused via a microscopic object 
lens 220 for optic-fiber input into a signal optical 
fiber 230. At the same time, the reference beam Ir is 
sent via a reflecting mirror 240 into a photo detector 
250 to generate heterodyne interference reference 
signals, which can be expressed as follows: 



in which F is the magnitude of the amplitude of the 
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heterodyne interference signal of the P wave component 
and the S wave component in the reference beam Ir- 

On the other hand, the signal beam Ig is outputted 
from the other end of the signal optical fiber 230, and 
5 diffuses into a multiple scattering medium 260 to form 
DPPDW. A detector fiber 270 similar to the signal 
optical fiber 23 0 is disposed in parallel at a separation 
distance (r) from the signal optical fiber 230 for 
detecting the DPPDW . In a photomult iplier tube (PMT) 280 , 

10 the heterodyne interference signals are generated, 
which are outputted through a signal amplifier 290 and 
a band-pass filter 300 to an electrical signal processor 
310 (e.g., a phase-lock amplifier) for amplitude and 
phase measurements. 

15 During the aforesaid process, when the polarized 

photon pair with mutually parallel polarization 
directions of the signal beam Ig enter into the multiple 
scattering medium 260, although the coherence and the 
degree of polarization (DOP) of most of the polarized 

20 photon pairs are decorrelated due toaseriesof collision 
events, during the process of detecting the optical 
interference signals , some of the polarized photon pairs 
can still maintain their correlation such that a 
meaningful heterodyne interference signal is generated, 

25 the main frequency of the interference electrical signal 

being Ao) = 0}p-0)^ . Therefore, after filtering through 
the band-pass filter 300 having a center frequency of 
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Aa> = cOp—co^ , other irrelevant noise signals can be clearly 
removed - 

Since only those of the polarized photon pairs which 
undergo relatively less collision events during passage 
5 through the multiple scattering medium can preserve 

relatively more polarization characteristics and 
directionality (or spatial coherence) for generating 
the heterodyne interference signals, the heterodyne 
interference technique according to this invention can 

10 be employed to screen diffused polarized photon pairs 
with a relatively low degree of scattering to form DPPDW . 
At the same time, they can also satisfy the diffusion 
equation. When the photon pairs with mutually parallel 
polarization directions are projected into a highly 

15 concentrated multiple scattering medium 260 and 
subsequently undergo photo detection, the intensity of 
the heterodyne interference signals thus generated can 
be expressed as follows: 

/^(Afu/) = Z)C + /cos(A£yr + AO) (2) 

20 in which ts^co^ (d^—o)^ is the beat frequency of the 

heterodyne interference signal; y is the amplitude of 

DPPDW; and AO is the phase delay. ia\ and //^ of the 

multiple scattering medium can be obtained from the phase 
delay and amplitude attenuation of DPPDW. Equation (2) 
25 can be further rewritten as: 
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/,(A^«>0=<Po — Rek^''""^'^} (3) 

r 



in which q>l is the rate of energy fluence of DPPDW; 

k^r is the real number portion of the wave number k2 of 

DPPDW, k2 =^2r +^"^2i ' ^ - V-1 ; and Atj^ corresponds to the 
absorption and scattering characteristics of the signal 
beam in the multiple scattering medium. 

^2r=[3//„fe+//J]^ (4) 

From Equation (3), different separation distances 

/"o and r can be obtained, and Ar = r-r(,. The amplitude 

attenuation of the signal beam in the multiple scattering 
medium is: 
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in which {Iq,I) are the intensities of the heterodyne 

interference signals at (r^.r) , respectively . Similarly, 
AO is the phase delay and can be expressed as follows: 
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in which k^i is the imaginary part of the wave number 



^2 of DPPDW, which corresponds to the scattering and 
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absorption characteristics of the multiple scattering 

medium. ( / /^a ) are, respectively, the reduced 

scattering coefficient and the absorption coefficient 
of DPPDW in the multiple scattering medium, n is the 
refraction index of the medium, and c is the speed of 
light. From Equations (4) and (7), the following 
equations are obtained: 

3w(Aa>) 



(9) 



2c 

10 Therefore, the characteristics ^ l^'s i l^a ^ of the 

multiple scattering medium can be obtained from the 
amplitude attenuation arid the phase delay of the 
heterodyne interference signals measured by the 
apparatus shown in Figure 1 . 

15 As shown in Figure 3, a two- frequency laser source 

400 and a detector 500 are displaced together relative 
to a multiple scattering medium 600 to perform 
two-dimensional (2-D) and three-dimensional {3-D) 

scanning. The spatial distribution of //^ and //^ is 

20 obtained from the amplitude and phase data of DPPDW. 

Then, the diffusion equation isusedto recover the ob j ect 
into an image. Details as to image recovery using the 
diffusion equation can be found in various literature, 
such as S. Fantini et al . , " Assessment of the size. 
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position, and optical properties of breast tumors in 
vivo by noninvasive optical methods," Appl . Opt, 37, 
1982-1989 (1998) ; and E.L. Heffer et al . , ^^Quantitative 
oximetry of breast tumors: a near- infrared method that 
5 identifies two optimal wavelengths for each tumor, " Appl . 

Opt. 41, 3827-3839 (2002) - The number of light sources 
and photo detectors can be adjusted suitably to obtain 
optimum imaging effects. 

Certainly, as apparent to those skilled to the art, 
10 the two- frequency coherent light source can be achieved 
using various constructions. For instance, a linear 
polarized single- frequency stabilized laser is used 
to emit a polarized beam, and the polarization angle 
thereof is adjusted through a polarization angle 

15 adjusting device , such as a A /2 wave plate - A polarized 

beam splitter splits the laser beam into orthogonal 
polarization directions, and both polarization laser 
beams are frequency shifted by passing through different 
frequency adjusting devices, such as acousto-optic 
20 modulators (AOM) . Hence, the frequencies of the two 
beams with mutually orthogonal polarization directions 

can have a discernible minor frequency difference A 

CO . Then, the beams are integrated once again to become 

a two- frequency coherent beam that includes the P wave 
25 and the S wave. The frequency adjusting device may be 

an electro-optic modulator or any other similar device . 
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Alternatively, by enabling a two- frequency laser with 
two different center wavelengths to undergo light 
synthesis, the polarized photon pair beam can also be 
projected into the multiple scattering medium by 
5 employing the aforesaid optical system so as to achieve 
a polarized photon pair for near- infrared spectroscopy 
(NIRS) . 

As shown in Figure 4, in another apparatus for 
implementing the method of this invention, rotation of 

X 

10 the azimuth angle through a — wave plate 110 can be used 

to adjust the azimuth angles of the P wave and S wave 
components of the two- frequency laser beam 100 and the 

X-Y coordinates so as to form an angle 0 (see Figure 5) - 

After splitting by a beam splitter 190, the signal beam 
15 Is is projected into the multiple scattering medium to 

be measured. Subsequently, during the process of 
detecting the optical interference signal of the signal 
beam Ig, the signal beam Ig is split once again by a 
polarized beam splitter 130 according to the 
20 polarization direction, and the mutually orthogonal 
linear polarized photon pair are respectively detected 
by photo detectors 140, 150 and signal processors 160, 
170, for conversion into heterodyne interference 
electrical signals, which can be expressed as follows: 

25 I X^o) t) = DC + y cos{ls.o} t + M>) (10) 
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Iy{l!s.(jot)^DC-ycos{^o)t + M>)'-- (11) 

By inputting synchronously the two heterodyne 
interference electrical signals obtained from 
conversion of the signal beam Ig that penetrated the 
5 multiple scattering medium to be measured into a 

differential amplifier 180, and by subtracting the two 
signals and amplifying the result for output, a balanced 
detector can be formed: 

A/ = /^ -ly =2ycos{^(Dt + M>) (12) 

10 The output of the differential amplifier 180 is fed 

to a signal processor 330. At the same time, the 
reference beam Ir is received by a polarizer 195 and 
a photo detector 320 from the beam splitter 190, and 
is subsequently inputted into the signal processor 330 

15 so as to obtain amplitude attenuation and phase delay 

of the heterodyne interference signals and to enhance 
the detection sensitivity and signal-to-noise ratio. 

Further, as the absorption coefficients //^ of 

hemoglobin in the blood in an oxygenated state (Hb02) 
20 and a deoxygenated state (Hb) are obviously different 

with respect to near infrared wavelengths ( >l, , Aj) (e.g., 

/I, =780nm ' X^^BSOnm) (see S. Fantini et al . , supra), 

changes in the hemoglobin saturation (Sa02) (see E.L. 
Heffer et al . , supra) in the blood can be obtained in 
25 real time from Equations (13) -(16). As relevant 
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theories can be found in the relevant literature 
mentioned hereinabove, the same will not be explained 
in detail herein. 

AMM=s^,AHb + s^,o,^0, (13) 

AMaM=^Hi>^b + e^,o,^b02 (14) 

(0-1.0^2 \ (15) 



✓2 rjuyj^, . 



^bO. -/> \ (16) 

K^HbO^^Hb ^HbO-i^Hb) 



in which iSHh ^ tsHbO^ ^ ts.fi^ are the variations of 

10 concentration of deoxyhemoglobin, oxyhemoglobin, and 
hemoglobin absorption coefficient , respectively. 

^Hb ^ ^Hb02 " " ^Hb02 extinction coefficients of 

hemoglobin to different wavelengths yl, and under 

deoxygenated and oxygenated states, respectively, and 
15 are known parameters. The amplitude attenuation and 
phase delay of DPPDW in the blood can thus be obtained 

therefrom. Then, //^ and ju^ of the blood are calculated 

according to Equations (5) - (9) . By utili zing different 

wavelengths and , the absorption coefficients 

20 Mai^) ^'^^ Mai^i) of hemoglobin can be obtained with 

precision, and the change in hemoglobin saturation in 
the blood can also be measured in real time. At the 
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same time, an image of the hemoglobin saturation in the 
blood can also be obtained by using spatial scanning 
in combination with the solution to the diffusion 
equation. This invention proposes a novel polarized 
5 photon pair near infrared spectroscopy, whereby 
measurements of the optical characteristics of a 
multiple scattering medium can be obtained . It is noted 
that results of experiments conducted-by the applicant 
are shown in Figures 6A and 6B, and it was found that 

10 the results match the anticipated results of Equations 
5 and 6 set forth herein [See Y.H. Chan et al . , ^^The 
measurement of optical properties of a multiple 
scattering medium based on diffused photon pair density 
wave". Advanced Biomedical and Clinical Diagnostic 

15 Systems, T. Vo-Dinh, W.S. Grundfest, D.A. Benaron, and 
G.E. Cohn, Proceedings of SPIE Vol. 4958,259-272 
(2003)] . In summary, by scanning the test object, in 
combination with the solution to the diffusion equation, 
the image of the object in a multiple scattering medium 

20 can be obtained. Moreover, by using a two- frequency 
laser source (preferably, two two- frequency lasers) 

having different center wavelengths, changes in Hb ^ 

HhO^ and Sa02, as well as the values of and IJ^J^^ 

that correspond to different wavelengths, can be 
25 measured with precision in real time. In case of a 
two- frequency , orthogonal linear polarized (or 
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circular polarized) photon pair, rather than a parallel 
polarized photon pair, propagating in the multiple 
scattering medium, by having a polarizer located between 
the scattering medium and the photomultiplier tube, a 
5 heterodyne signal can also be generated, and similar 
results can be achieved. While the present invention 
has been described in connection with what is considered 
the most practical and preferred embodiment, it is 
understood that this invention is not limited to the 
10 disclosed embodiment but is intended to cover various 
arrangements included within the spirit and scope of 
the broadest interpretation so as to encompass all such 
modifications and equivalent arrangements. 
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